Phonon lineshapes in the vortex state of the phonon-mediated superconductor 

YNi 2 B 2 C 
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We present an inelastic neutron scattering study of phonon lineshapes in the vortex state of the 
type-II superconductor YNi 2 B 2 C. In a previous study [Phys. Rev. Lett. 101, 237002 (2008)] it was 
shown that certain phonons exhibit a clear signature of the superconducting gap 2A on entering the 
superconducting state. Our interest was to find out whether or not the lineshape of such phonons 
reflects the inhomogeneous nature of the vortex state induced by a magnetic field smaller than the 
upper critical field B C 2 We found that this is indeed the case because the observed phonon lineshapes 
can be well described by a model considering the phonon as a local probe of the spatial variation 
of the superconducting gap. We found that even at B = 3T, where the inter-vortex distance is less 
than 300 A, the phonon lineshape still shows evidence for a variation of the gap. 

PACS numbers: 74.25.Kc, 78.70.Nx, 63.20.kd,74.25.Uv 



I. INTRODUCTION 

Many years ago, Axe and Shiranei were the first to 
show that phonons in phonon-mediated superconductors 
might exhibit marked changes of their frequency and line- 
shape if the compound is cooled below the superconduct- 
ing transition temperature T c . The effects are particu- 
larly pronounced if the phonon frequency is comparable 
to the energy of the superconducting gap 2A. For a long 
time, it was thought that the phonon lineshape can be 
well described by a Lorentzian above and below T c and 
the superconductivity- induced change of the phonon line- 
shape is just a change of the width of the Lorentzian. The 
first observation of a clearly non-Lorentzian lineshape be- 
low T c was made by Kawano et ali^ on the conventional 
superconductor YNi 2 B 2 C (T c = 15.2 K). The low tem- 
perature lineshape of certain phonons looked so unusual 
that the authors of this study thought to have observed 
the appearance of an additional excitation below T c . In 
a later study- 3 - a variety of strongly non-Lorentzian line- 
shapes were reported for the same compound. In particu- 
lar, it was shown that certain phonons exhibit a step-like 
increase of their intensity right at the superconducting 
gap energy 2A. It was shown that all the non-Lorentzian 
lineshapes are semi-quantitatively predicted by a theory 
proposed by Allen et alA 

In the present study, we address the question what 
happens to the phonon lineshape if a magnetic field is 
applied in a type II superconductor. It is quite obvi- 
ous to expect that the phonon lineshape will revert to 
that observed above T c if the magnetic field is so high as 
to completely suppress superconductivity, i.e. at fields 
larger than the upper critical field B C 2. If, on the other 
hand, the applied field is very small, i.e. less than B c \, 
the magnetic field does not penetrate into the supercon- 
ductor and hence no change is expected for the phonon 
properties. Therefore, we deal with fields B c \ < B < B C 2 



in the following where the field enters into the supercon- 
ductor in the form of vortices. In the core of the vortices, 
the superconducting gap is completely suppressed but 
gradually recovers its value without a magnetic field with 
increasing distance from the vortex core. This means 
that the vortex state of a type-II superconductor is inho- 
mogeneous with respect to the size of 2 A. In our exper- 
iment, we wanted to find out whether or not the spatial 
variation of the gap is reflected in the phonon lineshape. 

Theory is at present of little help to predict how 
phonons will respond to the inhomogeneity of the vor- 
tex state. In the most basic theory, i.e. in the harmonic 
approximation, phonons are represented by plane waves 
traveling through a perfectly periodic atomic lattice with- 
out damping. That is to say, not only is the phonon 
lifetime infinite, but there is also no length scale in the 
theory. An infinite lifetime is certainly unrealistic due 
to a variety of reasons, such as phonon-phonon coupling, 
i.e. anharmonicity, and in metals also electron-phonon 
coupling. This has led to the development of the quasi- 
harmonic approximation where phonons do have a finite 
lifetime, which can be seen in experiment as a broadening 
of the phonon lines. However, there is as yet no theory to 
deal with the length scale on which a phonon probes the 
microscopic properties of the solid. In some cases, one 
might get an estimate of the length scale by calculating 
the mean free path a; of a phonon from its lifetime and 
its group velocity via x = I/7 x v g , where 7 is the in- 
verse phonon lifetime and v g = dui / dq the group velocity. 
Both quantities can be obtained from inelastic neutron 
scattering measurements. For a strong coupling acous- 
tic phonon in Niobium, e.g., one obtains x = 700 A 5 , 
whereas for modes without strong coupling x is of the or- 
der of /im. This approach, however, is not applicable to 
zone boundary modes and many optic phonons because 
of their vanishing group velocity, and in particular not to 
the zone boundary mode studied in our experiment. 
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The mode chosen in our study was the endpoint of the 
transverse acoustic (TA) branch along the (110) direc- 
tion, i.e. at q = (0.5,0.5,0) (M-point). In a previous 
studj* 3 - we had found that the lineshape of this phonon 
shows a very steep intensity increase right at the energy 
of the superconducting gap 2A and therefore, a smearing 
or a shift of the gap could be easily detected in measure- 
ments in an applied field. Magnetic field dependent mea- 
surements have already been reported in Ref. |2j, using a 
TA phonon at q = (0.55,0,0) (see Sect. EVJ). The line- 
shape of this phonon does show pronounced changes on 
entering the superconducting state and therefore lends 
itself to a study of the influence of a magnetic field as 
well. However, the superconducting energy gap cannot 
as easily be extracted from the lineshape of this phonon 
as in the case of the M-point phonon 3 , and all the more 
any shift or smearing of the gap. 

In our experiment, we measured the evolution of the 
phonon lineshape of the M-point phonon as a function 
of applied magnetic fields 0T < B < 9T (B c2 = 9.3 T 
at T = 2K 6]). We show that at small fields, i.e. large 
inter-vortex distances, the phonon lineshape clearly re- 
flects the spatial variation of the superconducting gap in 
the crystal lattice. With increasing field strength, the 
superconductivity-induced deformation of the phonon 
lineshape is more and more reduced by the high den- 
sity of vortices, which makes it progressively difficult to 
obtain an unambiguous result. Our data obtained at 
B = 3 T, i.e. inter- vortex distances less than 300 A , still 
show evidence for a spatial variation of the gap, but in a 
less convincing way than at low fields. 



II. EXPERIMENTAL 

The neutron scattering experiments were performed on 
the cold triple-axis spectrometer 4F2 at the ORPHEE 
reactor at LLB, Saclay, (two runs, 1 and 2) and on the 
thermal triple-axis spectrometer PUMA at the research 
reactor FRM II in Munich. Double focusing pyrolytic 
graphite monochromators and analyzers were employed 
in both cases. At PUMA a fixed analyzer energy of 
14.7 meV allowed us to use a graphite filter in the scat- 
tered beam to suppress higher orders. Due to techni- 
cal limitations of the double focusing monochromator 
on 4F2, scans could be done using a fixed analyzer en- 
ergy of 14.7 me V only up to an energy transfer of 9meV, 
whereas for larger energy transfers the incident energy 
was kept constant at Ei = 23.7 meV. Therefore, the 
data had to be corrected for the k/ dependence of the 
analyzer efficiency. On PUMA, data were collected for 
fields 0.25 T < B < LOT, whereas data for higher fields 
were collected on the 4F2 spectrometer. The magnetic 
field was applied perpendicular to the (hhO) — (001) scat- 
tering plane, i.e. along the (110) direction. Changes in 
the applied magnetic field were always done in the non- 
superconducting state, i.e. at T > T c = 15.2 K (field 
cooled procedure). For each round of experiments zero 
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FIG. 1: (color online) Raw data obtained on the TAS 4F2 
of a transverse acoustic phonon at Q = (0.5,0.5,7) (run 1). 
Energy scans were measured at (a) B — and (b) 9T well 
below and above the superconducting transition temperature 
in zero field T c = 15.2 K. The vertical dashed line in (a) 
denotes the value of the superconducting gap at T = 2 K. 
Solid lines are fits to the data using a Lorentzian convoluted 
with the experimental resolution on an estimated background 
(dotted lines). The dashed line in (a) is a guide to the eye. 
The size of the shaded area in (a) was used to measure the 
strength of the lineshape deformation (see text and Fig. [4} 



field data were taken as well at base temperature and at 
a temperature slightly above the superconducting transi- 
tion temperature T c — 15.2 K. In order to achieve an ad- 
equate counting statistics, individual scans were repeated 
several times. Since it was clear that a spatial variation 
of the superconducting gap would show up at energies be- 
low ?s 8 meV, a large fraction of the total counting time 
for each scan was allotted to the low energy region. Un- 
avoidably, this strategy led to relatively large statistical 
errors at the upper end of the energy range investigated. 

The single crystal sample of YM2B2C weighing 2.26 g 
was grown by the floating zone method using the n B 
isotope to avoid strong neutron absorption^. The speci- 
men was mounted in a 10 T Oxford cryostat at 4F2 and a 
cryogen-free 7 T magnet on PUMA employing a standard 
closed-cycle refrigerator allowing measurement down to 
T = 2 K and 3 K, respectively. Measurements were car- 
ried out in the (hhO) — (001) scattering plane. The wave 
vectors are given in reciprocal lattice units (rlu) of (2n/a 
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2?r/& 2w/c), where a = b = 3.51 A and c = 10.53 A. 



III. RESULTS 

Energy scans performed on the TAS 4F2 in zero mag- 
netic field and B = 9T are shown in Fig. [TJ The 
zero field data taken at T = 2 K exhibit the abrupt in- 
crease of the phonon intensity at the superconducting 
gap energy 2 A = 5.7 meV as expected from our previous 
measurements^ albeit this feature is slightly washed out 
when compared to the results reported in Ref. H due to 
a slightly lower energy resolution of this instrument. On 
the other hand, the data taken on the TAS PUMA show 
again a very sharp feature at E = 2A, which allowed us 
to unambiguously detect the influence of magnetic fields 
as small as 0.25 T (Fig. [2]). 

We do not find any statistically significant difference 
between scans at T = 20 K and 5 K at B = 9 T (Fig. [lb). 
We note that the applied field of 9 T exceeds the upper 
critical at T = 5K 6 . This result docs not come as a sur- 
prise but just demonstrates that the distorted phonon 
lineshape is linked to the presence of the superconduct- 
ing phase and the corresponding energy gap. We now 
turn to the results depicted in Fig. [5J i.e. obtained at 
fields B < 1 T. The scans performed at B = 0.25 T 
or 0.5 T show the same step- like increase of the inten- 
sity at E — 5.7 me V as that performed in zero field al- 
beit the step height is significantly reduced. The line- 
shape in these scans can be very well reproduced by a 
superposition of the lineshapes observed in zero field at 
T = 3K and 17 K. In other words, the lineshape can be 
very well described by the assumption that the phonons 
sample volume fractions with zero or the full gap, re- 
spectively. The same model is fairly successful for the 
1 T-data as well, although a close look reveals that the 
step at E = 5.7 me V is somewhat smeared out. A more 
elaborate model to describe these data as well as those 
obtained at higher fields will be presented in the next sec- 
tion. At higher fields (Fig. a superposition of zero field 
data taken at T = 2K and 17 K, respectively, gives only 
an unsatisfactory account of the observed lineshape. This 
can be seen more easily comparing the superconductivity- 
induced changes, i.e. the difference between neutron 
counts in the normal and superconducting state, where 
each of the two data sets shown in Fig.(3t were subtracted 
from the counts measured at T — 17 K (Fig. [3Ji). Such 
a discrepancy had to be expected from the fact that in 
high fields, the gap will nowhere reach the full zero field 
value. 

In order to establish a quantitative relationship be- 
tween the strength of the superconductivity-induced line- 
shape deformation and the magnetic field, we evaluated 
the integrated difference of neutron counts between the 
lineshapes above T c and at very low T as shown in Fig.[T^ 
as gray shaded area. The results of this analysis are de- 
picted in Fig. 2J Somewhat surprisingly, the reduction of 
the lineshape deformation depends in a highly non-linear 
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FIG. 2: (color online) Raw data of energy scans of the trans- 
verse acoustic phonon at Q = (0.5, 0.5, 7) obtained on the 
TAS PUMA, (a) Comparison of zero field and B = 0.5 T 
data. Energy scans at (b) B = 0.25 T, (c) 0.5 T and (d) 1.0 T 
at T = 3K (dots). The open circles in (b)-(d) correspond to 
superpositions of the raw data at T — 3K and 17 K in zero 
field (see panel (a)) with the ratios given in the legend of each 
panel. 
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FIG. 3: (color online) Raw data of energy scans obtained 
on the TAS 4F2 of a transverse acoustic phonon at Q = 
(0.5,0.5,7) (run 2). Energy scans were measured in (a) zero 
field, (b) B = 2.0 T, (c) 3.0 T at T = 2K and T = 17 K 
(only zero field), (d) Enlargement of the low energy region 
of panel (c) including data taken at T — 12 K in zero field 
(open squares). The open circles in (b)-(c) correspond to su- 
perpositions of the raw data at T — 3 K and 17 K in zero field 
with the ratios given in the legend for each panel, (d) Dif- 
ference counts between normal and superconducting phase at 
small energy transfers, for which the data shown in (c) (same 
labeling) were subtracted from zero field data at T = 17 K 
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FIG. 4: (color online) Strength of the phonon lineshape defor- 
mation normalized to its maximal value at low temperature 
as function of magnetic field. The solid line is based on the 
vortex model as described in the text. Inset: The dashed line 
represents the magnetic field dependent lineshape deforma- 
tion of the acoustic phonon at q = (0.55, 0, 0) (reproduced 
from Fig. 4 of Ref. 0, see text). The solid line was calculated 
from our vortex model, but with a larger coherence length of 
^ = 125 A = v2 x 88 A corresponding to the reduction of the 
upper critical field by a factor of 2 from from B C 2 ~ 9.3 T at 
T = 2K to B c2 = 4.7T at T = 5.5 K. 



way on the magnetic field. Possible reasons for this ob- 
servation will be discussed in the next Section. 



IV. DISCUSSION 

We now want to discuss our results in the light of what 
is known for the vortex state of a type-II superconductor 
such as YM2B2C. Here, an applied magnetic field B en- 
ters the superconducting sample in the form of flux lines, 
so-called vortices, for T < T c and B c \ < B < B C 2- The 
superconducting gap in the center of a vortex is zero and 
the coherence length £ defines the radius around a vor- 
tex core where A reaches half of its maximum value. If 
a phonon would be a completely local probe, the exper- 
imental spectral weight distribution in the vortex phase 
should reflect the volume fractions having specific gap 
values. Starting from this assumption, we calculated the 
volume fractions with specific gap values using a value 
of £ = 88 A as reported in the literature^ and the func- 
tional form of the dependence of the gap on the distance 
from the vortex core as published in textbooks 9 . The 
magnetic field enters into this calculation via the inter- 
vortex distance a vortex , which is given by^£ 
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FIG. 5: (color online) Difference between neutron counts in 
the normal and superconducting state (dots) at B = 0.25 (a), 
0.5 (b), 1.0 (c) and 3.0 T (d) compared to predictions based 
on a vortex model (solid lines, see text). Experimental data 
were taken on the TAS's PUMA (a)(b)(c) and 4F2 (d). 



In the next step, we summarized the contributions from 
each volume fraction using the lineshapes observed for 
specific gap values in zero field in a previous experiment 3 . 
We note that these data, although collected with a dif- 
ferent set-up, can be very well transferred to the actual 
experiment after proper normalization (see Appendix). 

Results obtained in this way are shown in Fig. For 
the sake of an easier assessment of our approach, we 
did not plot the observed and calculated intensity distri- 
butions but instead the superconductivity-induced line- 
shape deformations for various field values. Apparently, 
our approach works quite well in the whole field range 
explored in our experiment. The good agreement even 
for B = 3 T is somewhat surprising in view of the fact 
that our model is oversimplified for short inter-vortex 
distances: it is known from the literature that the dis- 
tance dependence of the superconducting gap from the 
vortex core derived for isolated vortices is significantly 
modified for high vortex densities^. More specifically, 
it is expected that the gap size is more strongly reduced 
half-way to the next vortex core than calculated from 
the functional form for isolated vortices. Attempts to 
include such an effect into our model led to a serious 
disagreement with experiment. Later, we became aware 
that the good agreement between model and experiment 
at B = 3 T (and 2 T as well) is somewhat accidental 
because our model is somewhat unrealistic in another 
aspect: it assumes a field-independent value of the vor- 
tex diameter £ whereas investigations by muon spin res- 
onance have provided evidence for a drastic reduction of 
£ at high magnetic fields 12 . However, taking the value of 
£ = 40 A reported ini 2 . for B = 3 T at face value leads 
to a massive overestimation of the superconductivity- 
induced lineshape deformation: the effect should be as 
strong as calculated for £ = 88 A at B = 0.5 T (Fig. [3b), 
which is evidently incorrect. This discrepancy might be 
due to the fact the magnetic fields were applied along 
(001) in Ref . [l2l but along (110) in our work. We did not 
try to deduce a field dependence of £ from our data in 
order not to overstretch the information contained in our 
data, but concluded that there is insufficient knowledge 
of the properties of the vortex state in YNi 2 B 2 C at high 
fields for a realistic simulation of the phonon lineshape. 
An important qualitative conclusion can be drawn from 
our high field data nevertheless: the pronounced smear- 
ing of the feature related to the superconducting energy 
gap observed at B = 3 T indicates that the phonon line- 
shape reflects the inhomogeneous nature of the vortex 
state even for an inter- vortex distance as small as 300 A. 

A variation of the vortex diameter as function of the 
applied field from £ = 100 A for very small fields to 40 A 
at high fields seems to be unrealistic but on the other 
hand, our results do confirm such an effect qualitatively 
as evident from the highly non-linear influence of the 
magnetic field on the lineshape deformation. The mag- 
netic field dependence of the superconductivity-induced 
effect strength based on the structure of an isolated vor- 
tex is depicted as solid line in Fig. 2] Taking into account 
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that in regions with overlapping vortices the maximum 
gap value is further reduced^ the effect should fall off 
with increasing magnetic field even faster. However, our 
results show a larger effect for B > IT. This discrep- 
ancy can be qualitatively understood by a reduction of £ 
causing generally increased gap values compared to our 
model and a corresponding stronger deformation of the 
phonon lineshape. 

As mentioned in the introduction, the magnetic field 
dependence of the lineshape of another strong coupling 
phonon in YN12B2C was already investigated in Ref. 0- 
We pointed out that a shift or a smearing of the super- 
conducting gap cannot be easily derived from the line- 
shape of this phonon, because an evaluation of the su- 
perconducting gap relies heavily on the theory proposed 
by Allen et alj^, which, however, describes the observed 
effects only semi-quantitatively. A more summary in- 
formation can be obtained nevertheless: in the inset of 
Fig. [4] we compare the reported magnetic field dependent 
intensity of what the authors of Ref. |2j considered as a 
new peak with the strength of the lineshape deformation 
computed from our model outlined above. In order to 
account for the reduced value of the upper critical field 
at T = 5.5 K, the temperature at which the experiments 
reported in Ref. [U were performed, we used a value of 
£ = 125 A instead of 88 Ail A gain, we find a weaker 
reduction of the effect as function of magnetic field than 
expected for a constant value of £. Thus, the results for 
the TA phonon at q = (0.55, 0, 0) corroborate our conclu- 
sion that there it is indeed a decrease of £ as a function 
of magnetic field. 



100 



J2 50 
c 

O 

o 

o 
c 

CD 
CD 



-50 



40 



20 



c 

o 
o 

CD 

o 

CD 
I 



-20 



-40 



B = 0T, T = 3 K 
• PUMA'10 
O PUMA'05 

(devided by 4) i' 

» i 



J 1 T 1 , 1 j |_ 




(b) B=0T, T=2(3) K 

■ 4F2'11 
O PUMA'05 
(devided by 



4 

V Hi 



12 



energy (meV) 



FIG. 6: (color online) Difference between neutron counts in 
the normal and the superconducting state in zero magnetic 
field as observed in three different experiments. All data were 
properly normalized, (a) Data were taken on the PUMA 
spectrometer in the current investigation (filled symbols, with 
the sample mounted inside a superconducting magnet) and a 
previous experiment (open symbols, no magnet), (b) Data 
were taken in the early experiment on PUMA (open symbols) 
or the 4F2 spectrometer with the sample inside a lOT-magnet 
(full symbols). 



V. CONCLUSIONS 

Acknowledgments 



We have reported the magnetic field dependence of the 
lineshape of a certain phonon in YM2B2C, which exhibits 
a clear signature of the superconducting gap when the 
sample is cooled in zero field below T c . The observed line- 
shapes can be very well reproduced by a model based on 
the assumption that the phonon acts as a local probe of 
the inhomogeneous gap distribution in the vortex phase. 
The model is quite convincing for fields up to B < IT. 
Even at B = 3 T, where the inter-vortex distance is less 
than 300 A, the strong smearing of the superconductivity- 
induced feature indicates that the phonon resolves the 
gap values locally. A quantitative comparison between 
calculation and experiment at high fields is complicated 
by the fact that the vortex diameter seems to be field de- 
pendent in YN12B2C. Our analysis corroborates conclu- 
sions drawn from muon spin resonance dataii that the 
vortex diameter shrinks considerably in high magnetic 
fields in YNi 2 B 2 C. 
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VI. APPENDIX 

As explained in our discussion, the curves shown in 
Fig. [5] were calculated using observed phonon lineshapes 
for specific gap values. More precisely, we used zero field 
data taken at different temperatures whereby the sharp 
step in the intensity distribution allowed us to associate 
each temperature T <T C with a distinct gap value. How- 
ever, the limited beam time assigned to the experiments 
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in a magnetic field did not allow us to collect these data 
with the set-up used for the results reported in this pa- 
per. Therefore, we resorted to use data from a previ- 
ous experiment 3 performed without a magnet. We note 
that this experiment involved very long counting times 
and therefore led to good counting statistics. In Fig. [6] 
we demonstrate that the superconductivity-induced de- 
formation of the phonon lineshape observed in the ex- 
periments on the PUMA spectrometer is very much the 
same in different sample environments after proper nor- 



malization. The agreement between the data taken in 
2005 on the PUMA spectrometer and those in 2011 on 
the 4F2 spectrometer is not as good but still very sat- 
isfactory up to an energy transfer of about 8 meV. At 
larger energy transfers, there are systematic deviations 
between the PUMA and the 4F2 data, probably because 
the 4F2 instrument did not allow us to use a constant 
kf throughout the upper part of the scan. Therefore, 
our model calculations for the results obtained on 4F2 
spectrometer sit on firm ground only up to about 8 meV. 
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